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Abstract 

We  have  recently  measured  the  surface  phonon  spectrum  of  the  ordered  binary 
alloy  Cu3Au(001).  in  addition  to  the  Rayleigh  wave,  we  observe  a  higher  energy  mode 
which  may  be  interpreted  as  the  “folded"  Rayleigh  mode  or  "optical"  surface  mode.  Using 
a  harmonic  pair  potentials  fit  to  the  inelastic  neutron  scattering  data  of  Katar.o  et  a/,13 
we  find  that  the  force  constant  between  the  first  and  second  layer  Cu  atoms  must  be 
stiffened  by  20%  with  respect  to  the  bulk  value  in  order  to  match  the  "folded"  Rayleigh 
mode. 


Most  experimental  and  theoretical  surface  phonon  investiga¬ 
tions  to  date  have  focussed  on  semiconductor,  insulator  or  pure 
metal  surfaces.1  Few  experimental  phonon  studies  have  been  under¬ 
taken  on  the  surfaces  of  intermetallic  compounds  (ordering  binary 
alloys)  with  the  exception  of  two  recent  EELS  studies  on  bcc  NiAl 
(111),2  (110), 23  and  (100)2.  Intermetallic  compounds  are  of  inter¬ 
est  due  to  their  high  strength  characteristics  at  elevated  tempera¬ 
tures.4  These  remarkable  alloys  also  exhibit  increased  hardness  as 
compared  to  either  their  disordered  phase  or  their  single  metal 
constituents. 

In  this  study  we  examine  the  Cu3Au(001)  surface,  which  ter¬ 
minates  with  a  50%  Cu-  50%  Au  atomic  composition  in  a  c(2x2) 
surface  arrangement.5  In  this  orientation  the  material  consists  of 
alternating  layers  of  50%  Cu-50%  Au  and  100%  Cu.  This  A3B  inter¬ 
metallic  compound,  like  Ni3AI,  has  the  minority  species  occupying 
the  corners  of  a  face-centered  cube  with  (LI  2)  symmetry.  Cu3Au  is 
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a  good  model  system  for  investigating  how  interatomic  forces  differ 
between  the  bulk  and  surface  regions  uf  a  intermetallic  compounds.5 

CU3AU  is  particularly  interesting  because  the  bulk  undergoes  a 
first  order  phase  transition  (Tc=663  K)  to  form  a  substitutionally 
disordered  phase  at  a  temperature  well  below  its  melting  tempera¬ 
ture  (1226  K).5  Recent  theories6  and  experiments7  of  surface  order- 
disorder  transitions  indicate  that  they  may  occur  continuously,  but 
may  have  some  first  order  character  prior  to  the  bulk  transition.5 
Although  the  symmetry  of  the  lattice  is  reduced  when  going  from  the 
ordered  state  to  the  disordered  state,  the  material  retains  an  fee 
lattice — in  contrast  to  a  surface  melting  transition.8  This  behavior 
permits  investigation  of  the  surface  vibrational  properties  both 
above  and  below  the  transition  temperature.  We  report  here  the 
room  temperature  surface  phonon  dispersion  relations  of  the  ordered 
state  of  Cu3Au(001)  along  the  <100>  azimuth.9  Helium  atom  scatter¬ 
ing  is  well  suited  for  this  project  as  it  has  excellent  sensitivity  to 
the  periodic  arrangement  of  the  surface,  offering  superb  "contrast" 
between  the  topmost  and  underlying  layers.  It  also  serves  as  both  a 
structural  probe,  via  elastic  scattering,  and  a  dynamical  one,  via 
inelastic  scattering 

The  high  resolution  helium  scattering  apparatus  will  be 
described  in  a  future  paper.10  A  short  description  of  it  now  follows. 
The  machine  consists  of  a  helium  nozzle  beam  source,  an  ultra  high 
vacuum  scattering  chamber  and  a  rotatable  t.me-of-flight  mass 
spectrometer.  The  nozzle  source  can  produce  a  Av/v  velocity  spread 
in  the  incident  helium  beam  of  1%.  It  is  mounted  on  a  closed  cycle 
helium  refrigerator,  allowing  the  helium  beam  energy  to  be  tuned 
from  8-63  meV.  The  defining  aperatures  set  the  beam  angular 
spread  to  0.3°.  The  scattering  chamber  contains  standard  UHV 
preparation  (Ar+  sputtering  gun)  and  surface  analysis  tools  (AES. 

XPS,  and  LEED).  The  detector  consists  of  a  rotatable  quadrupole 
mass  spectrometer  mounted  on  optical  rails,  permitting  the  crystal- 
ionizer  flight  length  to  be  adjusted  to  either  55  or  100  cm.  I  he 
detector  can  independently  rotate  in  the  plane  of  scattering  a  total 
of  40  degrees,  allowing  scattering  data  to  be  collected  with  fixed 
incident  kinematic  conditions.  An  optical  encoder  mounted  on  the 
rotating  detector  and  interfaced  to  the  computer  keeps  track  of  the 
detector  angle  to  0.01°. 


15  20  25  30  35 


Final  Angle  -  3f 

Figure  1:  A  diffraction  scan  along  the  <100>  direction.  The  solid  line  is  a  fit  to  the  data 


The  Cu3Au(001)  sample,  previously  used  by  Graham,11  was  placed  in 
UHV  without  further  polishing.  The  surface  was  cleaned  by  cycles  of 
500  V  Ar+  bombardment  and  annealing  up  to  660  K.  The  crystal  was 
determined  to  be  clean  and  compositionally  ordered  when  the  60  eV 
Cu  and  the  70  eV  Au  Auger  peaks12  gave  a  Cu6o/(Cu6o+Au7o)  ratio  of 
0.52.  After  the  crystal  was  cleaned,  it  was  annealed  overnight  at 
approximately  630  K  in  order  to  enhance  the  growth  of  large 
domains.7  The  temperature  was  lowered  to  470  K,  1-2  hours  prior  to 
scattering.  Surfaces  prepared  in  this  manner  gave  sharp  full  and 
superlattice  LEED  diffraction  and  helium  diffraction  spots,  indica¬ 
tive  of  the  c(2x2)  surface  arrangement.  A  diffraction  scan  along  the 
<100>  direction,  showing  both  the  specular  and  superlattice  peaks 
(open  circles),  is  plotted  in  Figure  1,  along  with  a  fit  to  the  data 
(solid  line)  which  incorporates  the  instrument  function  of  the 
machine.  Upon  heating  the  crystal  to  663  K,  the  superlattice 
diffraction  spots  disappear,  while  the  full  order  diffraction  spots 
were  still  observed.  Occasionally  the  LEED  and  helium  diffraction 
spots  revealed  splitting  of  both  the  full  and  superlattice  peaks 
aligned  along  the  <100>  direction,  as  previously  reported  by  Me  Rae 
and  Malic.5 

After  the  crystal  was  aligned  along  the  <100>  direction  using 
the  superlattice  diffraction  feature,  a  series  of  time-of-flight  (TOF) 
measurements  were  taken  using  various  incident  energies  ranging 
from  16-32  meV  with  the  ionizer  positioned  55  cm  away  from  the 


Figure  2:  Energy  Transformed  TOF  spectra  taken  with  an  incident  energy  of  31.8  meV 
and  an  incident  scattering  angle  of  28-4°.  The  inset  shows  the  scan  curves  for  the  five 
runs. 


crystal  surface.  In  Figure  2  are  plotted  several  energy  transformed 
TOF  measurements  using  a  31.8  meV  (kj»7.78  A-1)  helium  atom  beam, 
with  an  incident  scattering  angle  of  28.4°  and  final  angles  ranging 
from  21.9°  to  25.9°.  Also  shown,  in  the  inset  of  Figure  2,  are  the 
scan  curves  showing  the  possible  energy  and  momentum  transfers 
for  a  given  detector  position.  In  all  TOF  spectra,  the  Rayleigh  wave 
is  clearly  discernible,  along  with  an  intense  diffuse  elastic  signal. 

In  some  of  the  spectra  a  nearly  dispersionless  higher  energy  mode  is 
also  observed,  essentially  a  "folded"  Rayleigh  wave,  present  because 
of  the  smaller  surface  Brillouin  zone  (SBZ)  of  the  ordered  c(2x2) 
surface.  All  of  the  spectral  features  plotted  in  Figure  3  are 


Surface  Phonon  Dispersion  Curves 
Cu3Au  <001  )-<l 00> 


Parallel  Momentum  Transfer  (A'1) 


Figure  3:  Dispersion  curves  for  Cu3Au(001)  along  r  -  M  '  in  the  extended  zone 
scheme.  Data  points  with  different  symbols  were  taken  on  different  days  and  correspond 
to  different  crystal  preparations.  Different  incident  energies  and  angles  were  used,  and 
the  crystal  was  approximately  room  temperature. 


displayed  in  the  extended  zone  scheme  along  the  r  -  M'  direction.  A 
total  of  about  50  TOF  runs  were  taken.  All  of  the  collected  points 
are  shown  in  a  folded  zone  scheme  in  Figure  4.  These  plots  also 
show  two  model  calculations  which  will  be  described  below. 

To  analyze  our  data  we  began  by  modeling  the  bulk  phonon  data. 
A  pair  potential  model  [similar  to  the  one  used  to  model  the  surface 
modes  on  NiAI(IIO)2]  was  used  to  fit  the  bulk  phonon  data  of  Katano 
et  a/.13  This  model  is  more  restrictive  than  a  Born-von  K&rm&n 
analysis,  and  ignors  noncentral  interactions.  The  two-body  poten¬ 
tials  consist  of  a  quadratic  part  [  -  fy"  x  (  |r|  -  |r0|  )2],  which  is 
essentially  the  force  constant,  and  a  linear  part  [  *  <t>j’  x  (  |r|  -  |r0|  )], 
which  is  a  tension  between  the  atoms.  Using  this  model,  a  nonlinear 
fit  to  the  bulk  phonon  data  for  CU3AU  was  performed,  incorporating 
second  nearest  neighbor  interactions.  The  results  of  this  seven 
parameter  fit  are  shown  in  Table  1. 
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*  indicates  the  modified  surface  tension  terms  used  in  Figure  4a. 
"indicates  the  modified  surface  force  constant  term  used  in  Figure  4b. 


To  characterize  the  surface  modes,  a  slab  calculation14  was 
performed  with  16  atomic  layers,  the  last  of  which  was  attached  to 
a  rigid  substrate.  Due  to  reflection  symmetry  through  the  chosen 
scattering  plane,  the  sagittal  modes  [shear  vertical  (SV)  and 
iongitudinal]  decouple  from  the  shear  horizontal  (SH)  modes. 
Selection  rules  dictate  that  our  experiment  was  sensitive  only  to 
sagittal  modes  Calculating  these  modes  required  diagonalizing  a 
64x64  dynamical  matrix  yielding  64  slab-normalized  eigenvectors, 
of  which  only  a  few  are  strongly  localized  to  the  surface  region. 
Helium  scattering  is  most  sensitive  to  SV  polarization,  so  that  only 
modes  whose  squared  SV  amplitude  in  the  top  layer  exceeded  0.5 
were  considered  likely  to  be  observed  in  our  experiment. 

When  the  theoretical  surface  derived  using  the  bulk  force 
constants  is  allowed  to  relax,15  ("driven"  by  the  nonzero  values  of 
<t>,'),  the  Cu  atoms  shift  upwards  0.2  A  out  of  the  surface  plane,  while 
the  Au  positions  are  unchanged.  Hence  the  surface  "rumples".  How¬ 
ever,  experimental  studies  of  the  atomic  geometry  of  Cu3Au(001)16 
and  0.5  ML  Au/Cu(001)17  seem  to  show  that  the  Au  is  displaced  out¬ 
ward  from  the  Cu  by  0.1  A,  differing  from  the  initial  theoretical 
prediction.  To  compensate  for  the  geometric  discrepancy  between 


this  simple  model’s  prediction  and  the  experimentally  observed 
rumpling,  the  magnitudes  of  pit,’  and  02b’  were  decreased  to  produce 
the  experimentally  observed  structure  (see  Table  1).  A  prediction 
for  the  surface  phonon  data  using  these  modified  tensions,  Ois'  and 
02s’.  is  shown  in  Figure  4a.  In  this  slab  calculation  the  Rayleigh 
mode  agrees  well  with  experiment,  but  the  upper  mode  is  too  low  in 
energy. 

To  fit  the  higher  energy  mode,  the  force  constant  02b"  between 
the  first  layer  Cu  atoms  and  the  second  layer  Cu  atoms  was  stiff¬ 
ened  by  20%.  (Figure  4b).  Stiffening  02S"  only  slightly  affects  the 
Rayleigh  wave,  while  it  greatly  influences  the  optical  mode.  This 
can  be  understood  by  considering  the  phonon  displacement  vectors. 
At  Q=  M‘,  the  Rayleigh  wave  displacement  vectors  involve  SV 
motion  of  the  first  layer  Au  atoms,  while  the  first  layer  Cu  atoms 
remain  essentially  stationary.  However,  for  the  optical  mode,  the 
first  and  second  layer  Cu  atoms  vibrate  normal  to  the  surface  plane 
180°  out  of  phase  with  each  other,  while  the  first  layer  Au  atoms 
vibrate  in  plane.  Hence,  the  optical  mode  at  Q=  M’  is  more  sensi¬ 
tive  to  force  constant  changes  involving  the  first  and  second  layer 
Cu  atoms. 


a) 
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Figure  4:  Lattice  dynamical  calculations  for  various  geometries  and  force  constant 
changes  (see  text),  a)  relaxed  positions  corresponding  to  the  experimental  data1 61 7 
and  using  the  bulk  force  constants.13  b)  same  relaxed  positions  as  in  (a)  but  with  a 
20%  stiffened  force  constant  connecting  first  and  second  layer  Cu  atoms. 
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The  20%  force  constant  stiffening  between  the  first  and 
second  layer  Cu  atoms  resembles  the  amount  of  stiffening  needed  to 
match  the  EELS  experimental  data  for  the  Rayleigh  wave  on 
Cu(001).18  In  fact,  the  force  constants  in  the  disordered  state  of 
CU3AU  are  almost  equal  to  the  force  constants  in  bulk  Cu.19  How¬ 
ever,  in  the  ordered  phase  of  bulk  CU3AU,  the  bulk  Au-Cu  first  near¬ 
est  neighbor  force  constant,  bib",  (which  is  about  3/2  stronger  than 
the  Au-Cu  bulk  force  constant  of  disordered  CU3AU),  is  five  times 
the  strength  of  the  bulk  Cu-Cu  force  constant,  <t>2t>"-  This  suggests 
that  the  surface  vibrational  properties  of  the  CusAufOOl)  surface  in 
the  disordered  state  may  be  very  similar  to  the  Cu(001)  surface. 

Of  further  interest  will  be  the  change  in  the  phonon  spectrum 
as  the  phase  transition  is  approached.  Recent  neutron  scattering 
work20  indicates  that  a  spinodal  transition  may  occik  below  Tc, 
(Ts=Tc-24  K),  since  the  bulk  phonon  modes  near  the  r  point  seem  to 
broaden  and  merge  into  one  peak,  starting  near  T=641  K.  Whether 
similar  behavior  is  seen  at  the  surface  will  be  explored  in  further 
experiments. 

In  conclusion,  we  have  presented  phonon  dispersion  results 
along  the  <100>  direction  for  the  ordered  Cu3Au(001)  surface,  where 
two  surface  modes  are  detected.  Using  a  central  potential  force 
constant  fitting  procedure  we  have  concluded  that  a  20%  stiffening 
of  the  force  constant  between  the  first  and  second  iuyer  Cu  atoms  is 
needed,  and  that  slight  relaxations  of  the  first  layer  nearest  neigh¬ 
bor  tensions  are  necessary  to  match  the  experimental  data.  Ab 
initio  calculations  would  be  useful  in  both  the  ordered  and  disor¬ 
dered  states  to  exDiain  the  change  in  bonding  as  this  phase  transi¬ 
tion  occurs. 
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